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Abstract. In this communication we introduce a pair of coupled continuum 
equations to model overlayer growth with evaporation-accretion due to thermal or 
mechanical agitations of the substrate. We gain insight into the dynamics of growth 
via one-loop perturbative techniques. This allows us to analyze our numerical data. 
We conclude that there is a crossover behaviour from a roughening regime to a very 
long-time, large length scale smoothening regime. 



PACS numbers: 71.20, 71.20c 
1. Introduction 

The dynamics of surface growth by atomic deposition have been the focus of interest 
over recent years (Q-f|). Several theoretical attempts (||) at the understanding of 
kinetic roughening have been made, through discrete and continuum models, motivated 
by experiments. Roughening is often an inevitable part of surface formation, so that an 
understanding of the surface morphologies has a crucial part to play in the many vital 
applications of this field. However, not much attention has been paid to the phenomenon 
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of smootfiening by thermal effects like evaporation. The physical picture is clear : the 
vibration of the substrate, of thermal or mechanical origin may smoothen the overlayer 
surface by transferring weakly bonded atoms on its bumps or mounds to available surface 
grooves. The picture is similar, but certainly not equivalent, to smoothening of granular 
surfaces by avalanches. This is in contrast to the surface diffusion term in continuum 
models, which arises because of the internal rearrangement of the bonded atoms in order 
to minimize the chemical potential. Avalanche smoothening has been recently studied 
in some detail using coupled continuum equations by Biswas et al ( [[10|] ). 

The notion that we shall borrow and adapt from the work of Biswas et al is that 
the dynamics of atoms in surface growth with evaporation-accretion is well described 
by the competition between the collective dynamics or relaxation of bonded atoms (in 
order to minimize the chemical potential) and the dynamics of free atoms diffusing on 
the surface. The bare surface of bonded, and therefore relatively immobile atoms, will 
be described by the local height h(r, t) above the substrate. Across this surface the gas 
of unbonded evaporated atoms diffuse until they are captured in an available groove. 
This gas of atoms will be characterized by its density p(r, t) just above the bare surface. 
A similar model has been discussed by Sanyal et al [1ITJ. However, the transfer term in 
the equations was rather different from ours and the possibility of smoothening was not 
addressed in that work. 

The usual practice for probing temporal or spatial roughness is to study the 
asymptotic behaviour of correlation functions like (h(r_,t)h(r!_,t')) in space or time via 
a single Fourier transform. Only one of the variables : space or time, is integrated over 
in Fourier space and the relevant scaling relations are used to determine the critical 
exponents that govern this behaviour. However, as pointed out by Biswas et al (||10||), 
this leads to ambiguities for those problems where there may be more than one scaling 
lengths. In such cases, the double Fourier transform provides a much deeper insight. To 
clarify this point, let us introduce the salient features of such a study : 

The connected two point self-correlation function of the local variable A(r, t), which can 
either be h(r, t) or p(r, t) is defined as : 



The scaling hypothesis implies that we have, in the absence of spatial anisotropy, 



S(r-r!,t-t') 



(A(r } t)A(r!,t>)) 



(A(r,t))(A(r f ,t')) 
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S(R, 0) 



R 
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R = r-r', R=\R 



and R 



oo 



for the saturated surface with t > t s , where t s is the saturation time, and 



S(0, r) 



T =\t 



r 



oo 



In general, 




for both R, r — > oo 
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The scaling function $ is universal, a is the roughness and z = a//3 the dynamical 
exponent. For the single Fourier transforms, 

S{k,T = 0) ~ AT 1 " 2 " for k^O 



and 



S{R 



0,u) 



10 



-1-2/3 



for 



(2) 



For the double Fourier transform, in case we assume strong scaling. That is, existence 
of single length and time scales, consequently a single dynamical exponent z. 

S{k, w) ~ uj- 1 ^ 1 - 2 * V for k, u -> 

which gives, 

S(k,u = 0) ~ jfe- 1 -**-* for fc^O 



and 

S(A; = 0,cj) ~ for (3) 

It is important to examine the single Fourier transforms in equations @. For the 
calculation of S(k, r = 0) we need to take the saturated surface after a long time, but 
for S(R = 0, oj) we need to take the entire growing surface but locally. Thus a is related 
to the saturated and (3 to the growing surface. However, the double Fourier transform 
requires information of both the growing and saturated surfaces, both locally and at 
large length scales. Biswas et al (||10||) correctly argue that, in case there are more 
than one length or time scales associated with the process, the double Fourier transform 
should provide a much clearer picture. 



2. The Statistical Model 

3. The statistical model for atomic deposition 

Atomic deposition has many features in common with granular deposition. The added 
feature is atomic binding. In the usual deposition geometry, a randomly fluctuating flux 
of atoms is incident on a substrate. Atoms deposit on the surface of the substrate and 
diffuse along it to minimize the energy. A cloud of unbonded atoms envelope this deposit 
and continuously exchange atoms with it through evaporation and re-deposition. 

While non-equilibrium growth has been extensively studied by coarse-grained 
classical stochastic equations ([0), it is not obvious a priori that the microscopic 
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energetic constraints relevant to atomic surfaces would automatically be satisfied by 
largely heuristic classical terms. In an earlier communication ( ||11]| ) we had presented 
electronic energy calculations in support of our model of surface growth. 

Among various physical processes which have been taken into account in models of 
growing interfaces, surface diffusion has been considered as the most important process 
involved. One such model involves the linear fourth-order Mullins- Herring continuum 
equation ([|l|, [13|]) supported by the discrete model of Wolf and Villain (WV) 



dh(r, t)/dt = -D V 4 /i(r, t) + r)(r, t) (4) 

where h(r,t) is the height of the interface from some mean height (h(r,t)) and f]{r,t) 
represents Gaussian white noise as usual. This equation yields a large roughness 
exponent a = 1.5 in d—1. 

In an earlier communication ([|ll|]) we had presented a model to look at the effect 
of desorption or evaporation on relatively immobile atoms which are bonded to the 
surface. A cloud of mobile atoms above the surface arise both from the impinging 
atomic beam and from evaporation caused by atoms knocked out of the surface by 
thermal or mechanical disturbances. These are described by their local density p(r,t). 
We propose a new class of growth equations with an explicit coupling between the profile 
of "bonded" atoms represented by the local height of the surface /i(x, t), and "mobile" 
atoms on the surface represented by their local density p(r,t). Our equations read: 



dh(r, t)/dt = - D h V 4 h(r, t) - T + Vh (r, t) 

dp(r, t)/dt = D p V 2 p(r, t) + T (5) 
where the transfer term T is given by : 



r = v 



V 2 h(r,t) 1 -Q(V 2 h(r,t)) - pp(r,t) V 2 h{r,t) ®{V 2 h{r } t)) (6) 



where O(x) is 1 for x > and otherwise. We describe in what follows the meaning of 
the above terms. 

(i) The fourth-order term in the equation(^j) describes surface diffusion of bonded 
atoms; this is the usual WV (|14||) term where represents a diffusivity. In the 



continuum picture surface diffusion is a manifestation of the collective motion of the 
bonded atoms leading to a shape rearrangement in order to minimize the chemical 
potential. The leading to this term is the gradient of the local chemical potential, 
which is assumed to be proportional to the local curvature. This assumption was 
shown to be valid provided we invoke the Locality Principle of Heine (||15||). 



(ii) The flowing atoms are neither bonded to one another nor to atoms on the surface. 
The first term in equation (|5|) describes normal, as opposed to surface diffusion of 
these mobile atoms, where the corresponding current is the gradient of the density. 
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(iii) The first term in the transfer term T, equation (0), describes spontaneous 
generation of mobile atoms on the surface through evaporation or clesorption. This 
could be due, for example, to thermal disturbances. We have assumed that it is 
easier to eject atoms weakly bonded at spiky mounds on the surface. For thermal 
ejection, for example, v is a measure of the substrate temperature. The Theta 
function ensures that it is easier to thermally eject atoms bonded on negative 
curvatures, i.e. on bumps or mounds on the surface. We have assumed that the 
rate of evaporation is proportional to the negative curvature on a mound. 

(iv) The second term in T, equation represents condensation, whereby mobile 
atoms accumulate and accrete preferentially at points of positive curvature, i.e. at 
the bottom of deep grooves where atomic coordination is large and hence bonding 
is strong. This term is obviously also proportional to the local density of mobile 
atoms. 

(v) Finally the last term is a Gaussian white noise characterized by its width 

(r ]h (r,t)r ]h (r!,t')) = A 2 h 5{r- r')8(t - t') 

We assume that growth occurs on a flat substrate; this and the absence of a preferred 
direction leads to the absence of anisotropy in space. In this model we have ignored the 
effects of Schowbel barriers. 

We can visualize the following sequence of processes: first, the mobile atoms diffuse 
(V 2 p) in the cloud above the surface. This is followed by the preferential conversion of 
these atoms into the bonded species at points of high positive curvature (p\ V 2 h\) on the 
surface such as mounds and grooves. The term v\ V 2 /i| models the effect of evaporation, 
leading to a dynamical exchange, at regions of high negative curvature, between bonded 
and unbonded atoms. However, the action of the V 4 /i term is to stabilize the formation 
of mounds and grooves, so ultimately the overwhelming effect is a competition between 
roughening and smoothening of the surface. Figure [I] illustrates the effect of the terms 
in our model. 



4. Numerical and theoretical analysis 

To start with we shall analyze a model in (1+1) dimension. The substrate is one- 
dimensional, so that the position of the overlayer on the substrate is located by a single 
position variable x. The coupled equations become : 



dh(x.t) _ d A h(x,t) d 2 h(x,t) , x d 2 h(xA) 

dp(x,t) d 2 p(x,t) d 2 h(x,t) d 2 h(x,t) 

~rjf~ = p ~cW~ + V ~te 2 ^ P[X,t) dx 2 (7) 
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Figure 1. A pictorial depiction of the model 



The Heaviside step functions in equation (§) introduces a complexity in our equations 
as far as analytic investigations are concerned. We shall follow the remarks of Biswas 
et al (|[10|) an d invoke a suitable representation of the Heaviside function as an infinite 
series. In that case, the equation (H) can be thought of as : 



_ d 2 h(x,t) , .d 2 h(x,t) 

T = » , 2 - , 2 + ... 

dx z dx z 

~ fd 2 h(x,t)\ n ~ . Jd 2 h{x,t)\ n 

We should note that the above expansion is not well-defined, as the coefficients of 
the expansion may themselves be large or divergent. As in fllU| , we shall still go 
forward in the spirit of self-consistency, i.e. subject to numerical verification. The 
Heaviside function introduces non-linearities. One way to gain some insight is to carry 
out a Hartree type mean-field approximation and replace the non-linearities by their 
expectation values. This leads to equations : 



dh(x, t) 
dt 



d A h(x,t) A d 2 h(x,t) , .d 2 h(x,t) 



dx 4 



dx 2 



dx 2 



dp(x, t) 
dt 



^ d 2 p(x,t) A d 2 h(x,t) A , .d 2 h(x,t) 

D n n : 1 + v — - ft p{x,t)- 



dx 1 



dx 2 



dx 2 



(9) 




where, /t = c/i and v = (l — c)v. Similar approximations have been studied by Bouchaud 
et al (0]). We expect that in some regime our equations (0) will reproduce the mean 
field results suggested by equations @. 

Figure |2] shows a part of the growing rough height profile for = D p — 1, fi — 1 and 
v = 0.01. The heights have been scaled in order to bring out the detailed features for 
comparison. We note that with increasing time short length-scale features slowly die 
out and mounds and grooves spanning longer lengths are formed. Figure |3] shows the 
variation of the root-mean square height deviation and density as a function of time. 
Both quantities show saturation with time. The full dynamical description then involves 
both the growing and the saturated profiles. 

The double Fourier transform is defined as : 

h(k,u) = jd X JdteM-Hkx-^h^t) 
The Green functions of the linear part of the equations are : 



8h(k, to) 
5i](k', uj') 



G h {k,uo) 5{k + k') 5{uo + uo') 
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Figure 3. The height (left) and density (right) root mean square deviations as 
functions of time, showing saturation 



/ 8p(k,u) 

\6r}(k',uj> 



G p (k,uj) 8(k + k') 5(u + u/) 



(10) 



The correlation functions are given by : 

Shikju) = (h(k,u) h(—k, —lj)) 
S P (k,uj) = (p(k,u) p(-k,-uj)) 



11) 



Now from equations (0), retaining only the linear terms, we obtain the two Green 
functions for h and p, 



9h(k,u 
9p(k,u 



) - (- 



-i uj + Dhk — vk A 
-l 



(-i uj + D p k'' 



(12) 



In order to describe the scaling behaviour of the deposition process, we usually 
define the following scaling indeces, from the correlation functions : 



S h (x -x',t-t') = (h{x,t)h(x',t')) - (h(x,t))(h(x',t')) 
S p (x-x',t-t') = (p(x,t)p(x',t')) - (p(x,t))(p(x',t')) 



S h (x, 0) ~ \x\ 2ah S p (x, 0) ~ \x\ 2a " as |x| -> oo 
S h (0,t) ~ \t\ 2l3h S p (0,t) ~ \t\ 2Ph as \t\ -> oo 



Study of a pair of coupled continuum equations modeling surface growth 
In general : 



9 



S h (x,t)~ \x\ 2a » F h {\t\/\x\*») 
S p (x,t)~ \x\ 2a » F p (\t\/\x\ z ») 

The scaling functions F h (£) and F p (£) are assumed to be universal and the indeces a^, a p 
and Zh = ah/Ph, z p = a p / (3 P are the roughness and dynamical exponents. Within the 
strong scaling hypothesis Zh = z p and there exists a single time scale and a distance 
scale. 



In the presence of two time scales, there is only a weak scaling [HJ hypothesis available 
to us : 



Here, Z/, 7^ z p . Absence of strong scaling means that the exponents and a p may 
become functions of k. 



4-1- Scaling analysis for the height-height correlations 



Let us look at the one loop diagram for the self-energy [|10 



E h (k,u) = fi 2 [ I ^- G h (k - q,u - cS) k 2 {k - q) 2 S p (q,u') 
2 f dq f duj 1 



2tt J 2vr 



1 



-i {u - u') + (k — qY - u(k - q) 2 + E h (k - q,w - u') 
k\k-q) 2 \ X p (q,u') 



(13) 



The integrand over the internal momentum q has a factor g (i+2a P ) ; which causes 



an infra-red divergence in the integral. As discussed by [10|, we take care of the infra-red 
divergence by introducing a lower cut-off ko <C 1. When we talk about "small values 
of momenta", we have to mean k ~ k Q . The argument for smaller values of k will have 
follow a different track. 

Note also that because of the term g - ( 1+2a p) in the integrand, the main contribution 
of the integral comes from small values of q. Now, for "small" internal momenta 
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q,w>-^o^--N q,w' 
k,w^ j>^> k,w k,w-->< 




k-q, w-w' 



E h (k,w) 



j>> k,w 

k-q, w-w' 

^k s w 

k,w ->A <>-> k,w 




E p (k,w) 



-> G h ( k ' w ) 

-> G p (k,w) 

if 2 

(xp V h vertex 



->0<- 



S h (k,w) 



S (k,w) 



>^-> vv h vertex 

Figure 4. One loop scattering diagrams for the self-energy for the Green functions. 
The scattering vertices are shown at the bottom of the figure 



k 3> q ~ k ; to' ~ g 2h , we can replace G h (k — q,cu — cu') ~ (« + A; 4 — z//c 2 + a;)) -1 . 
We now look at u; = 0, and noting that E(fc, 0) ~ A; 2 * 1 , we note immediately that if z ft < 4, 
then the inverse Green function is dominated by the self-energy : 

G^(M)~£,(M) 
Substituting this back in the equation for the self-energy : 

MMJ - £(lfc,0) Ap k 2vr J 2rr q*+*>>, 
This gives : 



Since, by definition S ft (A;, 0) ~ fc 2h , it follows immediately that z h = 2 
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k-q-Q,w-w'-W 



k-Q,w-W 




[1 (k,Q,k-Q ;w,w',w-w') 



k-Q.w- 

Figure 5. Vertex renormalization 



Here we have considered only the bare one- loop diagram. In general, the scattering 
terms leads to a renormalization of the fi vertex through the introduction of vertex 
functions fi(k, q,k — q) for uj = 0. The renormalization is illustrated in figure [5[ For 
q — > 0, we may assume fi{k, 0, k) ~ k x ^ . Putting these back into the equation for 
the self-energy, we note that Zh is renormalized to Zh + 8 where 5 = The exact 

numerical values for Zh may then differ from 2. 

The one-loop correction to the height-height correlation function is shown in figure 
(^). We may immediately write : 

q) A S h {k -q,u- (J) S p (q, u') 



S h (k,u) 



1 



C0 2 + \E h (k,u)\ 2 

Let us first examine the terms on the right-hand side : 



0U 2 +\Z h (k,Lj)\ 2 




(k ~ q) 4 
2tt (k-qY+ 2a h 1 



q 



l+2o p 



oo 2 + \E p (q,u')\ 2 J \(u-uj') 2 + \Z h (k-q,UJ-uj')\< 
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k-q,w-w' 



k-q,w-w' 



q,w 



k,w 



k-q,w-w' \ 




q,W 



k.w 



/ k-q,w-w' 



Figure 6. One-loop diagrams for the height-height and density-density correlation 
functions 



Going back to the expression for the self-energy, we see that : 
E h (k,u) : 

where 



^ 4 
-i uj + T k 2 



r 2 

1 o 



Substituting this in the expression for the correlation. 

S h (k,u) = 



f r 4L8 

u 2 + L ° k 



T Q k 2 



u 2 + T 2 k\ 



to 2 + r^ 4 



(14) 



This is our main result from which various limits may be obtained. For instance, putting 
uj = : 



S h (k,u = 0) 



vlk 1 V r„P, 

~ k~ % for small k 



The numerical results are shown in figure [7]. The best fit slope was found to be 6.75. 
Now, using the expression 

S h (k,to = 0) ~ k~ {1+2ah+Zh) 

we obtain an equation : 



1 + 2 a h + z h = 6.75 



(15) 
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Figure 7. The Log-log plot of the double Fourier transformation S^kjU) — 0) vs k. 
For small k the slope is —1 — 2a/, — Zh — -6.75 

Again, putting k=0 and referring to figure [8] we get : 

Sh(0,ou) ~ uj~ 2 

Numerically this has been found to be 1.996. Using the expression, 
S h (k = 0,tu)~cu-( 1+2l3h+1 ^ 





ln(w) 



Figure 8. The Log-log plot of the double Fourier transformation Sh(k = 0,uj) vs u. 
For small u> the slope is —1 — 2/3/, — l/z/, ~ -1.996 



Study of a pair of coupled continuum equations modeling surface growth 14 
we get another equation, 

1 + 2/3 h + l/z h = 1.996 (16) 

Using the equations ([15]) and ([[£]) and Zh = cnh/Ph we estimate : 

a h = 1.2 f3 h = 0.36 and z h = 3.38 (17) 

Note that numerical estimate of Zh is greater than the one- loop estimate of 2. As 
discussed earlier, the one-loop estimate is perturbative and may be considerably modified 
by vertex renormalization, as well as higher order diagrams. 

We get further information from the full double Fourier transform. For example, 
for sufficiently small k, but w ^ we get, 

S h (k, u) ~ A{lo) + B{lo) k 2 




Figure 9. The Log-log plot of the double Fourier transformation Sh(k,oj) vs (top) ui 
for different values of k and (bottom) vs k for different values of u . 
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I I i I i I i I i I i I 

-5 -4 -3 -2 -1 

ln(w) 

Figure 10. The Log- log plot of the single Fourier transformation (top) Sh{k,t = 0) 
vs k, fitted to the slope —1 — 2 or/, = —3.746 and (bottom) Sh(x — 0, u) vs uj fitted 
to the slope -1 - 2f3 h = -1 .763 

As k — > this leads to a flattening of the curve to a u-dependent constant A(u). The 
value of k at which the flattening occurs goes on decreasing as uj decreases. Similarly, 
for / 0, for sufficiently small u> : 

S h (k,u) ~ C{k) +D{k) u 2 

As before, the value of uj at which the flattening occurs decreases as we decrease k. 
This crossover behaviour is seen in our numerical results. It is easy to observe that 
we would not have gleaned this information from single Fourier transforms, a point 
emphasized earlier in the work of Biswas et al (|nj). We may interpret this as a long- 
time, smoothening of the growing surface and is a characteristic of the evaporation- 
accretion process. 

In the range of k prior to flattening, we may derive : 
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S h (k,t = 0) = J^S h (k,u) 

~ J_ f dU + n 2 k 3-2a h+Zh f fa 

~ 2tt J w 2 + k**H ^ J {uj 2 + k 2 ^) 2 

~ A k' Zh + B k 3 ~ 2ah ~ 2zh 

Estimates of a>h and (3h are 1.37 and 0.38 respectively, consequently Zh = 3.61. Note 
that there is no reason why the estimates from single Fourier transforms and the double 
Fourier transform should agree exactly. For the single Fourier transforms we need 
information only about the saturated surface for Sh(k,t = 0) and the growing surface 
for Sh{x = 0,uj). For the double Fourier transform, we need information for both the 
saturated and growing surfaces for the same function. In case there are multiple length 
scales in the problem that information will be reflected in the double Fourier transform 

m 



From figure [It] in the k <C k region, flattening suggests that : 

Sh(k, t — 0) ~ const ~ k° 

Since Sh(k,t = 0) ~ k~( l+2ah \ the numerical results suggest that in this regime 
at = —0.5 and Zh = 0. In this regime, we have Sh(k,tu = 0) ~ fc-^+^h+Zh) ^ 
This crossover to a flattened regime for very small k is also seen in the numerical results 
of figure 0. 

4-2. Scaling relations for density- density correlations 

Again referring back to figure (§|), we can write an expression for the self-energy for the 
density Green function : 



E p (k, u) = /x 2 / p- I ^ G p {k -q,u- J) S h (q, u;')k 2 (k - qf + u 2 k A S h (k, u) 



.2 



dq f du' ( 1 \ 1 



2tt J 2tt \-i (iv - u') + EJk - q,u - u') I q 1+a h 



1 



\oj' 2 + \E h (q,uj')\ 2 J 

We may now carry out the integral over uj' and use the fact that the integrand over q has 
a factor q~( 1+a h) which ensures that only small q values contribute to the integral. This, 
combined with a lower cut-off for q, ko allows us to evaluate the q integral approximately. 
The procedure is almost exactly like the case of height-height correlations : 
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Figure 11. Log-log plot of the double Fourier transform S„(k,U! = 0) vs k. The best 
fit slope is : —1 — 2a p — z p = -6.49. 



E p (k,uj) 



T 2 k A 



-i uj + Tik 2 



For small k, we know that H p (k,uj = 0) ~ k Zp , and from above equation : 

E p {k,uj = 0) ~ A k 2 + B k 3 - 2ah ~ Zh 

With ah = 1-5 and Zh = 2, we obtain z p = 2 

For the density-density correlation function we get (again from figure ( 



S p (k,u) 



u 2 + \Y, p {k,u)\< 



v 2 k^ 

+ 2 T t.2z S h(k,U>) 
U) A + K p 



The first integral over q has a integrand with a factor g~( 1+a '') ) which makes sure that 
the main contribution to the integral comes only from small values of q. So in the 
integrand, we can replace k — q by q and carry out the u' integral, to obtain : 



S p (k,uj) ~ A yu 2 + 

Given zu = 2, z p = 2 and an = 1.5. 

S p (k,u = 0) ~ A;" 6 



-i 



to 2 + T 2 k* 



Tpk 2 
u 2 + k 
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Figure 12. Log-log plot of the double Fourier transform S p (k — 0,u>) vs to. The best 
fit slope is : -1 - 2f3 p -l/z p = -2.48. 



From figure [11], numerically we get an index of 6.49. 

In the cross-over regime when Zh = , ah = —0.5 and z p = 2: 



S p (k, uj — 0) 



fc(2z h +2z p ) 



k 



This is seen numerically as a flattening in the very low k regime in figure |TI 
Again, integrating over the variable uj we get : 



S p (k,t = 0) 



(TS/2) 



k z p(k Zh + k z p) 



Numerically we find the index to be 3.868. 

From the expression for S p {k,uj) we note that : 

S p (k = 0, uS) ~ w~ 2 

Numerically we obtain an index of 2.488. We may also carry out the integral over k to 
get : 
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The numerical prediction for the index is 2.05. The accompanying table summarizes 
our results. 
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1 t~\ /H /~\ 

illQcX 

(Single Loop) 


\ 1 1 t v i ot* 1 1 * ' i 1 
IN LLlllt/J. IL-cLl 

lIlQeX 




A tici l^r'l'ir'ci 1 
JrYLlaiy Llt^cLl 

lllUtiX 

(Single Loop) 


\f n TYi pri r* q 1 
IN Lllllt/i IL-dl 

1 v~» /H "x ^ 

iiiuex 


Tih(k, uj = 0) v s k 


2.00 


3.74 


Sp(fe, = 0) v s k 


2.00 


2.87 


Sh(k, t = 0) v s k 


4.00 


3.75 


S p (k, t = 0) vs k 


4.00 


3.89 


Sh(k = 0, u) vs to 


2.00 


1.996 


S p (k = 0, to) vs u 


2.00 


2.48 


Sh(k, uj = 0) v s k 


6.00 


6.65 


S p (k, uj = 0) v s k 


6.00 


6.49 



5. Conclusion 

We have studied the scaling behaviour of a set of coupled continuum equations describing 
surface growth in the presence of evaporation-accretion both numerically and within a 
one-loop perturbative approach. We notice a crossover from a roughening regime to 
smoothening at very large time, large length scale regimes. The one-loop estimates 
gives us an insight into our numerical results. We also notice that in order to study 
problems with more than one length-time scales, i.e. if we have weak scaling, it is 
essential to study the double Fourier transforms which probe both the growing and the 
saturated profiles together. This is in confirmation of the ideas set forward earlier by 
Biswas et al (PA). 
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